The Number of Same-Sex Marriages in a Perfectly Bisexual Population is Asymptotically Normal 
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Theorem: Consider a population of 2n men and 2n women where every individual is equally 
attracted to either sex and chooses his or her mate according to other criteria. Also assume 
that everyone gets married. Then the expectation of the random variable "Number of same-sex 
marriages" is 
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and its variance is 
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Furthermore, this random variable is asymptotically normal. 

Semi-Rigorous Proof: The probability generating function is (why?) 
^ ! 
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^ I By repeatedly applying the operation f{x) — )• xf'{x), plugging in x = 1, for n = 1 to n = 200, 

and fitting the data by a rational function of n, one gets in turn the mean, variance, and higher 
moments, from which one easily gets the moments about the mean and the normalized moments (a 
coefficients). 



^ I It turned out that the mean and variance are indeed as stated by the theorem. As for the higher 

moments, we will show that the first 14 normalized moments tend (as n — >• oo) to those of the 
normal distribution, and the readers can go on as far as they please. 
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http: //www. math. rutgers . edu/ "zeilberg/tokhniot /SameSexMarriages . 

More (computer-generated!) detailed output can be gotten from the front of this article 

http: //www. math. rutgers .edu/"/mamarim/mamarimhtml/ssm.html . 
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The normalized third moment (about the mean) is 
1/2 V2 
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the asymptotics of its square is: 
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Note that it goes to 0, as it should, since the odd moments of the normal distribution are all 0. 

The normalized fourth moment about the mean (alias kurtosis) is 

3 - 100 n + 650 - 1896 + 2632 n'^ - 1664 + 384 
' n2 (35 - 188 n + 348 - 256 + 64 n^) 

Its asymptotics is: 
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Note that it converges to 4!/(2^ • 2!) = 3, as it should, this being the fourth moment of the normal 
distribution. 

The normalized fifth moment about the mean can be viewed in 

http://www.math.rutgers.edu/"zeilberg/tokhniot/oSameSexMarriagesl, and the asymptoticS of its Square 

25 , 315 3501 7 37067 o 383273 „ ^, 

y " + ^ " + 12^ " + ^ " + W " + ^ ) • 
Note that it goes to 0, as it should, since the odd moments of the normal distribution are all 0. 

The normalized sixth moment about the mean can be viewed in 

http://www.math.rutgers.edu/~zeilberg/tokhiiiot/oSameSexMarriagesl . ItS asymptoticS is: 
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Note that it converges to 6!/(2^ • 3!) = 15, as it should, this being the sixth moment of the normal 
distribution. 

Etc. etc. (Sec http: //www. math. rutgers.edu/'zeilberg/tokhniot/oSameSexMarriagesl for evidence Up tO 

the 14*'^ moment, and you are welcome to modify http://www.math.rutgers.edu/"zeilberg/tokhniot/inSameSexMarriagesl 

and run it on your computer, if you want more evidence.) □ 
Comments 

1. It should be routine, once we found the average and variance as above (that even humans should 
be able to do easily), to derive a local limit law, by using Stirling's asymptotic formula. However, for 
some reason. Maple refuses to take the appropriate limit. (It works for the binomial distribution). 
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2. Another automatic, this time fully rigorous proof, can also be gotten using an extension of the 

methods in [Zl], using the Maple package 

http : / / www . math . rutgers . edu/~zeilberg/tokhniot/CLT. However, since the polynomials Pn (.x) 
are not closed-form, but rather satisfy a second-order linear recurrence (easily derived via the Zeil- 
berger algorithm) one should be able to do it (but CLT has to be slightly extended). 

3. The present approach follows [Z2] and the Maple package accompanying it: 
http : // www . math . rutgers . edu/ "zeilberg/ tokhniot /HISTABRUT. 

4. Not directly related to the present article, but nevertheless highly recommended, are references 
[AZl], [Z3], [Z4], and [Z5]. 

References 

[AZl] Moa Apagodu and Doron Zeilberger, Multi-variable Zeilberger and Almkvist-Zeilberger algo- 
rithms and the sharpening of Wilf-Zeilberger theory, Advances in Applied Mathematics 37 (2006) 
139-152. 

[Zl] Doron Zeilberger, The Automatic Central Limit Theorems Generator (and Much More!), "Ad- 
vances in Combinatorial Mathematics: Proceedings of the Waterloo Workshop in Computer Algebra 
2008 in honor of Georgy P. Egorychev", chapter 8, pp. 165-174, (I.Kotsireas, E.Zima, eds. Springer 
Verlag, 2009.) Available on-line from 

http : / / www . math . rutgers . edu/ zeilberg/mamarim/mamarimhtml/ georgy . html . 

[Z2] Doron Zeilberger, HISTABRUT: A Maple Package for Symbol- Crunching in Probability theory. 
The Personal Journal of Shalosh B. Ekhad and Doron Zeilberger, Aug. 25, 2010. (Journal Main 
page: http://www.math.rutgers.edu/ zeilberg/pj .html, direct url of article: 
http : //www. math. rutgers . edu/ zeilberg/mamarim/mamarimhtml/histabrut .html) 

[Z3] Doron Zeilberger, Andre's reflection proof generalized to the many-candidate ballot problem. 
Discrete Math 44 (1983), 325-326. 

[Z4] Doron Zeilberger, Kathy O'Hara's constructive proof of the unimodality of the Gaussian poly- 
nomials, Amer. Math. Monthly 96(1989), 590-602. 

[Z5] Doron Zeilberger, The Umbral Transfer- Matrix Method. I. Foundations Journal of Combina- 
torial Theory, Series A, 91 (2000), 451-463. 



3 



